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We study the vacuum, plane- wave Bianchi Vllh spacetimes described by the Lukash met- 
ric. Combining covariant with orthonormal frame techniques, we describe these models in 
terms of their irreducible kinematical and geometrical quantities. This covariant description 
is used to study analytically the response of the Lukash spacetime to linear perturbations. 
We find that the stability of the vacuum solution depends crucially on the background shear 
anisotropy. The stronger the deviation from the Hubble expansion, the more likely the overall 
linear instability of the model. Our analysis addresses rotational, shear and Weyl curvature 
perturbations and identifies conditions sufficient for the linear growth of these distortions. 
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1 Introduction 



There has long been an interest in the study of the spatially homogeneous Bianchi spacetimes 
and their cosmological applications to our understanding of singularities and of the observed level 
of isotropy in the universe. These studies analyse the problems within the manageable domain 
of ordinary differential equations and provide only a finite number of alternative cosmologies 
(see [1] and references therein). The most general Bianchi universes that contain the open 
Friedmann model as a special subcase are those of type Vllh- The late-time asymptotes for 
the non-tilted type Vllh spacetimes, with h ^ and a matter content that obeys the strong 
energy condition, evolve towards the vacuum plane-wave solution found by Doroshkevich et al 
and Lukash [2, 5] that is known as the Lukash metric. These spacetimes describe the most 
general effects of spatially homogeneous perturbations on open Friedmann universes; see for 
example [6]- [10]. The Lukash metric plays a guiding role in these investigations because of the 
subtle stability properties of isotropic expansion at late times in open universes. When the 
strong energy condition is obeyed, then isotropic expansion was found to be stable but not 
asymptotically stable at late times [7]- [10]. 

Traditionally, Bianchi spacetimes have been studied qualitatively, primarily by means of dy- 
namical system methods [11]-[16]. The same techniques also facilitate the analysis of the less 
well understood tilted Bianchi models, namely those where the fluid 4- velocity is no longer or- 
thogonal to the hypersurfaces of constant time [17]- [20]. In this paper we attempt an analytical 
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approach. In particular, we combine covariant and orthonormal frame techniques to provide a 
description of the Lukash spacetime of Bianchi type Vllh in terms of its irreducible kinematical 
and geometrical quantities. We then use the zero-order results to study inhomogeneous pertur- 
bations around the vacuum solution and discuss its linear stability. In so doing, we allow for the 
presence of a low-density matter component with a pressure-free equation of state. Our main 
interest is the evolution of perturbations in the kinematics and the geometry of the vacuum 
model and its linear response to these distortions. We find that stability depends primarily on 
the amount of the background shear anisotropy. In particular, our results show that the higher 
the background shear the more likely is the linear instability of the model. Another key factor is 
the relative orientations of the various kinematical and geometrical quantities. The positioning 
of the vorticity vector with respect to the principal axes of shear, for example, or the relative 
orientation between the shear and the spatial curvature eigenframes can also influence the linear 
stability of the Lukash spacetime. 

We also consider rotational perturbations and study their evolution relative to the back- 
ground volume expansion. Our analysis provides a condition for the linear growth of vortical 
distortions that depends primarily on the background shear anisotropy. We identify, in partic- 
ular, the minimum amount of shear necessary for the linear instability of the Lukash universe 
against vorticity perturbations. The level of background shear can also determine whether kine- 
matical anisotropies will remain bounded or dominate the linear expansion of the perturbed 
Lukash model. Additionally, by allowing for low-density dust matter, we find that its presence 
has no effect on the vorticity and shear of the perturbed spacetime. In contrast, the introduction 
of a material component can have an effect on the linear Weyl anisotropy of the model. More 
specifically, the presence of even a low-density dust fluid seems enough to ensure that the Lukash 
solution will diverge from its original plane- wave nature at late times. 

2 The Lukash plane-wave attractor 

The Bianchi Vllh models belong to the non-exceptional family of the Behr class B spatially 
homogeneous spacetimes. The plane- wave Lukash solution is the late-time attractor of the 
Bianchi Vllh models for a broad range of initial date and matter properties. These vacuum 
spacetimes correspond to equilibrium points of the associated autonomous dynamical system 
and are self-similar [21]-[24]. The line element of the Lukash metric takes the form 

ds 2 = -dt 2 + t 2 dx 2 + t 2r e 2rx [(Ady + Bdzf + {Cdy + Adz) 2 ] , (1) 

where r is an arbitrary constant parameter in the range < r < 1, A = cost;, B = / _1 sini;, 
C = —fsinv and v = k(x + hit) [25, 26]. Note that / and k are constants related to r by 

k 2 

— (1 - / 2 ) 2 = 4r(l -r) and r 2 = hk 2 , (2) 

where h is the associated group parameter. As we shall see next, constraint (2) is the Lukash 
analogue of the Friedmann equation. We also point out that for r = 1 and f 2 = 1 the Lukash 
metric reduces to that of the empty Milne universe. 

3 Covariant description 

Consider a family of observers, with worldlines tangent to the timelike velocity field u a (nor- 
malised so that u a u a = — 1). The latter, together with the associated projection tensor h a b = 
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9ab + u a Ub, introduces a local 1+3 threading of the spacetime into time and space. One can then 
decompose the various kinematical, dynamical and geometrical quantities into their respective 
irreducible parts and obtain a completely covariant description of the spacetime [27, 28]. 

3.1 Covariant variables 

The covariant formalism uses the irreducible kinematic quantities, the energy density and pres- 
sure of the matter fields and the gravito-electromagnetic tensors, instead of the metric which 
in itself does not provide a covariant description. The key equations are the Ricci and Bianchi 
identities, applied to the observers' 4-velocity, while Einstein's equations are incorporated via 
algebraic relations between the Ricci and the matter energy-momentum tensors. Thus, in the 
absence of matter and rotation, the plane-wave attractors of the Bianchi Vllh spacetimes are 
covariantly characterised by 1 

H = = p = q a = ir ab and u a = = u a , (3) 

while 

Q,a ab , E ab , H ab ^0. (4) 

Note that fi, p, q a and ir ab are respectively the energy density, the isotropic pressure, the heat 
flux and the anisotropic stresses of the matter, 0, a ab , uo a and ii a are the volume expansion, the 
shear, the vorticity and the acceleration, while E ab and H ab are the electric and magnetic parts 
of the Weyl tensor (C a bcd)- The latter have equal magnitudes and are orthogonal to each other, 
in accord with the Petrov type N nature of the Lukash solution. In other words, 

E 2 = H 2 and E ab H ab = , (5) 

where E 2 = E ab E ab /2 and H 2 = H ab H ab /2. The former of these constraints implies that 
CabcdC abcd = 0. The latter ensures that C abcd *C abcd = 0, where *C abcd = r] abqs C* s cd /2 and 7] abcd 
is the 4-dimensional alternating tensor. Note that the Weyl curvature invariant C abcd C abcd has 
been suggested and used as a measure of the gravitational entropy by several authors [29]-[33], 
but cannot on its own capture deviations from isotropic expansion in plane-wave spacetimes. 

3.2 Covariant equations 

The average volume expansion of the Lukash universe is described by the following version of 
the Raychaudhuri equation 

G = -±6 2 -2ct 2 , (6) 

where a 2 = a ab a ab /2 is the magnitude of the shear tensor. As usual, the expansion scalar is 
used to define an average scale factor (a) via the standard relation G/3 = a/ a. 

The absence of matter means that the Lukash spacetime is Ricci flat. The curvature of the 
spatial sections, however, is not zero. In particular, zero rotation ensures that the 3- Ricci tensor 
(Hab) is completely determined by its scalar and its symmetric and trace-free parts. These are 
given respectively by 

U = -|6 2 + 2a 2 , (7) 

Sab = -\®&ab + Vc{a° C b) + E ab , (8) 

throughout this article we employ a Lorentzian metric with signature ( — , +, +, +) and use geometrised units 
with c = 1 = 8irG. Consequently, all geometrical variables have physical dimensions that are integer powers of 
length. Also, Latin indices take the values 0,1,2,3 and Greek ones run from 1 to 3. 
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where S ab = TZ( ab ) = Tt( a b) ~ TZh ab /3. 2 As we will see below, the scalar 1Z is negative, which 
means that the model is spatially open. Note that expression (7) is the generalised Friedmann 
equation. 

In covariant terms, gravitational waves are described by the electric and the magnetic parts 
of the Weyl tensor. The latter obey a set of three coupled propagation equations, which are 
accompanied by an equal number of constraints. In the case of the Lukash plane-wave spacetime 
the evolution equations take the form 

E ab = -eE ab + cur\H ab + 3a c{a E c b) , (9) 
H ab = -@H ab -curlE ab + 3a c{a H c b) , (10) 

&ab = -§©0"<ib - E ab - (J c (a(J C b ). (11) 

The constraints, on the other hand, are 

B b a ab = |D a e, (12) 
B b E ab = e abc a b d H cd , (13) 
B b H ab = -e abc a b d E cd , (14) 

where e abc = r\ aba }U d is the spatial alternating tensor. In addition, the shear and the magnetic 
component of the Weyl tensor are directly related by 

Hab = c\n\a ab , (15) 

with cui\H ab = e cd ^ a D c H d b ^ by definition. Clearly, exactly analogous expressions define cvLi\E ab 
and cui\a ab . Note that the presence of (standing) gravitational waves is guaranteed by the 
non-zero values of both cm\E ab and cux\H ab . This is possible, despite the spatial homogeneity 
of the Lukash spacetime, because of the non-zero 3-curvature of the model. In other words, 
D c E ab , D c H ab ^ due to non-zero Christoffel symbols. 

4 Orthonormal-frame description 

The orthonormal frame formalism is an 1+3 decomposition of the EFE into evolution and 
constraint equations relative to the timelike vector field eo of an orthonormal frame {e a } [34]- 
[37]. In cosmological studies eo is the fundamental 4- velocity field, usually identified with the 
motion of the cosmic medium. In models that accept an isometry group, eo can also be chosen 
as the normal to the spacelike group orbits. 

4.1 Structure constants 

The Lukash solution belongs to the non-exceptional family of the Bianchi class B spacetimes. 
For these models the structure constants n a p and a a (with n a pa a = 0) take the form 3 

n a/3 = diag (0, n 2 , n 3 ) and a a = (a±, 0, 0) . (16) 

2 Angled brackets denote the symmetric and trace-free part of orthogonally projected tensors and the orthog- 
onally projected components of vectors. 

3 In the orthonormal formalism the spacetime metric is g a t = r/ a t = diag (— 1, 1, 1, 1) and the spatial frame 
vectors are {e a }. Greek indices are raised and lowered by means of the spatial metric g a p = S a p. 
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Moreover, the self-similarity of the Bianchi Vllh plane-wave attractor guarantees that the three 
non-zero components of the structure constants are given by [25] 

r k . kf 

ai = --, n 2 = — and n 3 = — , (17) 

where r, k and / are constants related with each other and with the group parameter of the 
model (see Eq. (2)). 

4.2 Kinematics 

Non-exceptional, non-tilted Bianchi class B spacetimes, like the Bianchi Vllh cosmologies, have 
a 12 = = (T13. Hence, the self similarity of the vacuum plane- wave attractor of these models 
means that the remaining components of the shear tensor are 

2(1-0 I-*" , fc(l"/ 2 ) 
011 = ^ . 0-22 = 0-33 = ^— and CJ23 = ^— . (18) 

On using the above, one finds that the magnitude of the shear tensor associated with the Lukash 
solution is 

„2_i n gp_ (l~r)(l + 2r) 

where a a p are the non-zero orthonormal frame components of a a b. On the other hand, the mean 
Hubble volume expansion of the Lukash universe is determined by the scalar 

1 + 2r 

Q = ^ r L - (20) 

The above ensures that the average scale factor obeys the simple power law a oc 

t (l+2r)/3 with 

< r < 1. As expected, for r — > 1 the scale factor evolution reduces to that of the Milne 
universe (i.e. a oc t). Note that the effect of the anisotropy is to reduce the expansion rate below 
that of the isotropic case. 

When measuring the average anisotropy of the expansion, it helps to introduce the following 
dimensionless and expansion-normalised shear parameter 

In the Lukash spacetime the scalars a 2 and O are given by (19) and (20) respectively. Using 
these expressions we obtain 

1 — r 

E = , (22) 

1 + 2r ' y 1 

This result reflects the model's self similarity, which guarantees that all the expansion-normalised 
dimensionless variables remain constant in time. Given that S > and < r < 1, we immedi- 
ately deduce that < S < 1, in accord with 1Z < in (7). Thus, although the shear anisotropy 
is not asymptotically stable (in the Lyapunov sense), it is stable in the sense that any deviations 
from isotropy never diverge [8]- [10]. Note that for r — > 1 the ^-parameter approaches zero and 
the expansion becomes isotropic. Maximum shear anisotropy (i.e. S — > 1), on the other hand, 
corresponds to the r = limit. 
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Definition (21) also provides an alternative expression for the Raychaudhuri equation of the 
Lukash solution. In particular, combining (6) and (21) one arrives at 

9 = -±6 2 (l + 2£) , (23) 

while the power-law evolution of the average scale factor now reads a oc 

t i/(i+2E)_ Xhug) inthe 

absence of any shear anisotropy we have a oc t, as in the Milne universe. For maximum shear 
anisotropy, on the other hand, we arrive at the familiar scale-factor evolution of the Kasner 
vacuum solutions (i.e. a oc i 1 / 3 ). Note that the deceleration parameter of the Lukash model is 
q = 2S, which means that q = when £ = and takes value q = 2 of the Kasner points as 
E -» 1. 

4.3 Spatial curvature 

In a vacuum, plane-wave Bianchi Vllh spacetime, the trace of the 3-Ricci tensor 7Z a p associated 
with the surfaces of constant time is 

n = -i(n 2 - n 3 ) 2 - 6a? , (24) 

where a±, ri2 and n% are given in Eq. (17). According to the above 1Z < always, which 
guarantees the hyperbolic geometry of the spatial sections. Note that we may use (17) to recast 
Eq. (24) as 

fc 2 (l-/ 2 ) 2 6r 2 _ 2r(l + 2r) 

Then, combining results (19), (20) and (24), one can show that that expression (7) (i.e. the 
Lukash analogue of the Friedmann equation) reduces to the constraint (2a). 

Spatial curvature anisotropics are described via the symmetric and trace- free tensor S a p. In 
the Lukash model, the only non-zero components of S a p are 

4r(l - r) _ fc 2 (l-/ 2 )(2 + / 2 ) fc 2 (l - / 2 )(1 + 2/ 2 ) 

611 _ §tj — ' 622 37^2 ' 633 y 2b ) 



523 = -^V^- ^) 



and 

£r(l_- f) 
ft 2 

According to (25)-(27), the spatial curvature of the model vanishes at the maximum shear limit, 
namely as r — > 0. At the other end, as r approaches unity, only the isotropic part of TZ a p 
survives. Recall that k 2 (l — f 2 ) = as r — > or 1 (see Eq. (2a)). 

4.4 Weyl curvature 

The only non-zero components of the Weyl curvature tensor associated with a vacuum, plane- 
wave Bianchi type Vllh spacetime are 

E22 = — -E33 = —#23 = — ^2f 2 J ~ 

and 

H - H -F _ fc(W 2 )(l~2r) 

H22 — — -n 3 3 — ^23 — • V 2y J 
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This means that the Weyl curvature is minimised near the limits r = 1 and r = 0, where 
k 2 (l — f 2 ) = 0. Also, when r = 1/2, which corresponds to £ = 1/4 (see Eq. (22)), we have 
H22 = = i?33 = E23 and the Weyl tensor has only one independent component. 

These relations between the Weyl tensor components also guarantee that the electric and 
magnetic Weyl tensors have equal magnitudes, although they are orthogonal to each other. More 
specifically, expressions (28) and (29) imply that 



E = H 



k 2 (l~f 
Apt 4 



2\2 



^(l + / 2 ) 2 + (l + 2r) 2 



(30) 



and that 

E a/3 H^ = 0. (31) 

As with the expansion anisotropy, it helps to measure the anisotropy of the Weyl field by 
means of the following expansion-normalised, dimensionless scalars 



W+ = 



W+ 

04" 



and 



W- = 



04" 



(32) 



where W± = E 2 ± H 2 by definition. The self-similarity of the Lukash solution guarantees that 
W+ is time independent, while the plane-wave nature of the model ensures that W- = (see 
constraint (30)). 



5 The perturbed Lukash solution 

The Lukash solution is the late-time attractor of the Bianchi Vllh spacetimes [8, 9], which are 
known to contain the open FRW universe as a special subcase. In this respect, studying the 
behaviour of the perturbed Lukash model could provide useful clues to the final stages of ever- 
expanding FRW cosmologies with /j + 3p > 0. If /J, + 3p < then the expansion will approach the 
FRW (for power-law inflationary behaviour [38]) or de Sitter universe in accord with the cosmic 
no-hair theorems. Similar effects can arise from the effects of higher-order curvature corrections 
to the Einstein-Hilbert lagrangian of general relativity [39]. 



5.1 Nonlinear equations 

Consider a perturbed vacuum Bianchi Vllh spacetime and allow for a low-density, pressure-free 
matter component. The nonlinear evolution of the latter is governed by the standard energy- 
density conservation law 

/i = -6/x. (33) 

When the matter component is in the form of dust, there is no acceleration and the only 
additional kinematic contribution comes from possible rotational disturbances. This means 
that the cosmological velocity field, which is identified with the motion of the matter, remains 
geodesic although it is allowed to rotate. Rotation is monitored by the propagation equation of 
the vorticity vector 

Co a = -|Bo; a + a ah uj b , (34) 

which also satisfies the constraint 

B a oj a = . (35) 
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According to (34), in addition to the expansion effect, which always reduces vorticity, there is a 
contribution due to the shear anisotropy. Note that for pressure-free matter there are no sources 
of rotation and vorticity remains zero if it was zero initially. 

In the presence of a pressureless fluid and vorticity the nonlinear Friedmann and Raychaud- 
huri equations respectively give 4 

U= -§0 2 (l-X-ft) -2lu 2 (36) 

and 

e = -ie 2 (l + 2£ + ift)+2u; 2 , (37) 

where O, = 3/u/0 2 is the density parameter. Similarly, the introduction of matter and rotation 
modifies the rest of the propagation formulae given in section 1.2 as follows 

E a b = ~QE ab + cur\H ab + 3a c ( a E b f - \na ab - u c e cd ( a E b ) d , (38) 
H ab = -@H ab -cuTlE ab + 3a c{a H b) d -u; c e cd{a H d } , (39) 

<?ab = -^QcTab- E ab - O c ( a V C b) ~ ^{a^b) , (40) 

while the associated constraints become 

H ab = curler^ + D (a u; b ) , (41) 

D b a ab = |D a 9 + curW a , (42) 

V b E ab = e abc a b d H cd -3H ab u b + lD af i, (43) 

B b H ab = -e abc a b d E cd + ?,E ab uj b + ^ a . (44) 

Note that vorticity affects the expansion and the shear evolution only at the nonlinear level, 
while it has a linear contribution in Eqs. (38), (39) and (41)-(44). 

5.2 Linear vortices 

The presence of matter means that one can identify the cosmological velocity field with that 
of the material component, which in turn gives physical substance to the idea of rotation. We 
measure the relative strength of rotational perturbations by means of the expansion-normalised 
dimensionless scalar 

G7=|, (45) 

with uj = (ojauj -) 1 / 2 . Taking the time derivative of to and using (23) and (34) we obtain the 
linear expression 

w = -\Q (l - 2S - 3£ afe nV) w , (46) 

where the tildas indicate background quantities. Here, T, ab = a ab /@ by definition and n a is the 
unit vector along the rotation axis (i.e. uj a = un a ). Written in an orthonormal frame the above 
reads 

d t w = -±6 (l - 2X) w + eX afS n a nPm , (47) 

4 In a low density perturbed model with Q -C 1, the linear Friedmann and Raychaudhuri equations retain their 
background functional form. The difference is that E is generally not constant, which means that the average 
scale factor of the perturbed Lukash model no longer obeys the simple power-law evolution given in section 3.2. 
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where S a/ g and n a , rip are the non-zero orthonormal frame components of T,^ and n a , 
respectively. The first term on the right-hand side of this equation describes the average evolution 
of w, while the second conveys the directional effects. The former increases (or decreases) w 
depending on whether £ is greater (or less) than 1/2. 5 Overall, linear vortices grow, relative 
to the average background expansion, when the following condition holds 

2S + 3± a pn a nP > 1 . (48) 

This condition implies that the growth of linear vortices also depends on the relative orientation 
between the background shear eigenframe and the rotation axis. Assuming that rotation takes 
place along the ei axis of the background orthonormal frame, we may use expressions (18a) 
and (22) to verify that condition (48) holds as long as r < 1/2. 6 In this case linear vortices 
grow relative to the background expansion as w oc t l ~ 2r . This means that the growth rate 
of w takes the maximum value w oc t at the r = limit, namely for maximum background 
shear. Alternatively, one may assume that the rotation axis lies along e2 or e%. Then, a 
similar calculation shows that linear vortices can never grow relative to the average background 
expansion. At best, w remains constant (when r — > 0). 

Following expressions (46) and (47) , the average shear distortion always increases the residual 
amount of rotation. Moreover, when condition (48) is fulfilled, the overall effect of the shear 
(including the direction dependent component £ Q| g) will also boost vorticity perturbations. In 
other words, as far as rotation is concerned, shear distortions can mimic the effects of matter 
pressure. Recall that in the presence of pressure vorticity does not necessarily decay with time 
(e.g. see [27, 40, 41]). Instead, for matter with a stiff enough equation of state rotation will 
increase despite the universal expansion. In our case, non-zero pressure with p = p(fi) means 
that Eq. (46) takes the form 



1 - 2S - 3 (c* + £ a6 nV) w , (49) 



where c 2 = dp/dfi is the square of the adiabatic sound speed. This demonstrates clearly the 
analogy between the shear and the pressure effects on rotation. For example, when E = 1 and 
the direction-dependent term on the right-hand side of the above is negligible, the shear effect 
on w is indistinguishable from that of a matter component with pj \x = 2/3. 



5.3 Linear shear anisotropies 

In a perturbed Bianchi Vllh model with low-density dust, Raychaudhuri's formula (see Eq. (37)) 
ensures that the linear expansion proceeds unaffected by the presence of matter or by rotational 
distortions. Also, expression (40) guarantees that, to linear order, vortical perturbations do not 
affect the evolution of the expansion-normalised shear parameter. Thus, ignoring rotational and 
matter effects, we take the time derivative of (21) and then use Eqs. (6), (8) and (11) to arrive 

5 The value E = 1/2, which corresponds to r = 1/4 (see expression (22)), indicates a Lukash-type spacetime 
with half the allowed amount of shear anisotropy. In the case of rotational distortions, E = 1/2 also indicates 
the point where the background expansion rate drops faster than the vorticity. Recall the decreasing effect of the 
shear on the average expansion scalar (see Eqs. (6), (23)). Therefore, at E = 1/2 the scalar vj starts to increase 
on average (i.e. excluding the direction dependent effect of the last term in (47)). 

6 The condition r < 1/2 on the metric parameter translates into the constraint E > 1/4 on the background 
shear anisotropy. Accordingly, only three quarters of the allowed Lukash backgrounds are potentially unstable to 
linear rotational perturbations. 
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at 

£ = -§6(l-£)£- ^S ab a ab . (50) 

The first term in the right-hand side of this equation describes the average linear evolution of 
£. Conversely, the last term of Eq. (50) describes directional effects and depends on the relative 
orientations of the S a b and a a b eigenframes. Relative to an orthonormal coordinate system the 
above reads 

d t £ = -f6(l-£)£- J^/*' 3 , (51) 

with S a/ 3 and o a p representing the non-zero orthonormal frame components of S a b and a a b 
respectively. 

Suppose that £ — > to zero order. This state corresponds to a background of minimum shear 
and 3-Ricci anisotropy (i.e. r — > 1) with & a p = = S a p, where tildas indicate the zero-order 
quantities. At this limit, which corresponds to a perturbed Milne universe, Eq. (50) takes the 
linear form 

d t £ = -§€>£. (52) 

As a result, £ oc a~ 4 oc t~ 4 to first order. In other words, when £ — > any linear expansion 
anisotropies that may occur will quickly disperse. 

When £ — > 1 we have r — > and the background has maximum shear anisotropy and zero 
3-Ricci curvature. This means that S a g = but & a p ^ (see expressions (18), (26) and (27)). 
Here it helps to introduce the auxiliary linear variable S = 1 — £ and rewrite Eq. (52) as follows 

dtS = |OS + ^-d aP S^ , (53) 

given that S — > 0. Note that whenever S grows the £-parameter decreases and vice-versa. Also 
note that although the background S remains bounded within the open interval (0, 1), this is 
not necessarily the case at the linear level. Indeed, using the zero-order relations (18) and the 
trace-free nature of S a p, we find that the last term on the right-hand side of (53) equals SSut 
in the r — ► limit. This in turn allows us to recast Eq. (53) as 

d t S=-^S + 3S 11 t, (54) 

where Sn is a component of the perturbed spatial Ricci tensor. The above means that a decrease 
in S and therefore an increase in the linear shear anisotropy when £ — > 1 is possible in principle. 
Suppose, for example, that the perturbed Sn component retains its background form, namely 
that Sn = -4r(l - r)/3t 2 with < r < 1. Then, Eq. (54) solves to give 



ft \ 4/3 
* = +3.(1 



to J 



4/3' 



(55) 



where Sq = S(to). Recalling that S = 1 — £ by definition, this leads to the following expression 
for the perturbed shear parameter: 

, 3r 2 (l-2r) / t\ 4/3 , , 

2 = i _ 3 r(l - r) + } n > - , (56) 



l + 2r \toJ 
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assuming that S = (1 — r)/(l + 2r). The latter means that the initial relation between the 
perturbed S and r has the background functional form. Given that r < 1, expression (56) 
implies that the perturbed ^-parameter can break through the S = 1 barrier provided that 

^ 4/ % (l-r)(l + 2r) 

to J > r(l-2r) • (57) 

For example, if r ~ 10 -2 , this will happen when t > 10 3 / 2 io- Note, however, that the smaller 
the value of r the longer it takes for the perturbed £ to cross through unity. 

It should be emphasised that S > 1 to first order implies that the linear 3-Ricci curvature 
of the perturbed Lukash model becomes positive (see Eq. (36) with Q = = u> a ). This is 
possible for maximum background shear (i.e. as r — ► 0), because the zero-order spatial Ricci 
tensor vanishes at that limit. So, in principle, the perturbed spacetime can have spatial sections 
with slightly positive curvature. Clearly, in this case the linear Lukash model perturbs away 
from the family of the Bianchi Vllh cosmologies, a fact demonstrated by the unbounded shear 
parameter. On the other hand, if we demand that the linear 3-curvature is never positive, then 
E will always remain bounded by unity. 



5.4 Linear Weyl anisotropy 

Consider the expansion-normalised dimensionless variable W_ defined in (32a). This scalar 
vanishes in the Lukash plane- wave background, which implies that a linear growth for W- is 
a sign of instability at that perturbative order. In addition, W_ directly determines the Weyl 
curvature invariant C abc d,C abcd . In this respect, the linear evolution of W_ also monitors the 
gravitational entropy of the perturbed Lukash universe and, to a certain extent, that of low- 
density open FRW models. 7 

For weakly rotating Lukash universes with a low-density dust component (i.e. when zu, <C 
1), the time derivative of definition (32b) leads to 

VV_ = -|(l-4E)eW- + ^ i (^ a6 curlF a5 + F a6 curl J B ab ) 

+ J^ca (E c b E ah - H c b H ab ) - ^ c e cda (E b d E ab - H b d H ab ) 

~^4^E ab , (58) 

on using Eqs. (37)- (39). The last two terms in the right-hand side of the above describe the effects 
of vorticity and matter respectively. Note that although matter does not directly contribute to 
the Weyl field, the later is not entirely arbitrary because of the contracted Bianchi identities. 
These are in a sense the field equations for the Weyl curvature and, among others, convey the 
matter effects on the propagation of the Weyl components (e.g. see [42]). On introducing an 

7 A complete description of the Weyl anisotropy of the perturbed vacuum Bianchi Vllh spacetime, also requires 
to study the linear evolution of the expansion normalised scalar W+ (see definition (32a)). Unlike W-, however, 
W+ has nonzero background value. This complicates further the linear study of W+ and allows only for relatively 
trivial analytic solutions of the associated propagation equation. It is conceivable that an improved version of the 
formalism presented here will be able to address the full Weyl anisotropy of the perturbed Lukash solution. 



11 



orthonormal frame, expression (58) reads 



d t W- = -l{l-AT l )@W- + ^ va (E^d^H v p + H^d^E v p) 



e 4 

3 

04"^"^ y- " p 1 " " py q4' 



- 0T (V^ " ^/^) " ^4 • ( 59 ) 

Again, we notice that the first term on the right-hand side describes the average evolution of 
W_, while the rest describe direction-dependent effects. The former depends crucially on the 
background shear anisotropy and it is reversed in sign at S = 1/4. This threshold corresponds 
to r = 1/2 and indicates the point where the background expansion rate starts decreasing 
faster than the average W-. Recall that at S = 1/4 the background Weyl field has only one 
independent component (see section 3.4). 

To proceed further we consider homogeneous perturbations in the Weyl field, namely that 
d^Eap = = d^Hap. This condition also monitors the large-scale behaviour of >V_ in the 
presence of inhomogeneities. In addition, we assume that E^E^ = E a pE^ v and HapH^y = 
H a pH^ u to linear order. 8 Then, the combined linear contribution of the second, third, fourth and 
fifth terms in the right-hand side of (59) is zero, while the sixth term reduces to —2(1 — r)W-/t. 
Also, employing the background relations between the Weyl tensor components given in (28) 
and (29), one can immediately show that the vorticity term in Eq. (59) vanishes to first order. 
Finally, on using the zero-order expressions (2), (18), (20), (28) and given that fi cx 

r (l+2r) for 

dust (see Eq. (33)), we find that 

» * pa/3 _ I 1 * p „ r(l-r)(l-2r) 2r 

W UapE ~ e 4 " 23 (TT2r)4 ' (60) 

to first order. Note that, according to this, the matter effect in Eq. (59) vanishes when either 
(T23 or ^23 is zero. On these grounds the linear expression (59) reduces to 

W_ = -2r^W_ + r(1 ( 7 + )( 2 1 r ^4 2r) ^- 2r . (61) 

where the parameter r varies in the open interval (0, 1) and C is a constant. This result implies 
that the expansion-normalised scalar W_ remains unchanged to linear order as r — > 0, that is 
for maximum background shear. The reason is that in the r = limit the Weyl anisotropy of 
the Lukash solution disappears (see Eqs. (28), (29)). The matter effect also vanishes near the 
minimum background shear limit (i.e. as r — > 1) and at the r = 1/2 threshold. The former of 
these two results is not surprising, since the Lukash solution decays to the Milne universe at the 
r = 1 limit. When r = 1/2, however, the matter effect is zero because E23 = at that point 
(see Eq. (29)). Recall that at the r = 1/2 threshold the Weyl tensor has only one independent 
component. 

When r -> we find that c?tW_ = 0, which implies that any deviations in W_ that may 
occur will remain constant. On the other hand, as r — ► 1 or at the r = 1/2 threshold the 

8 This assumption, which allows us to obtain analytic solutions for the linear evolution of W- , implies that the 
£11 and Hn components of the perturbed model vanish. This restriction already holds in the background. 
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expansion-normalised Weyl parameter decays as 

W- oc r 2 and W_ oc r 1 , (62) 

respectively. In general, ignoring the effect of matter leads to W_ oc t~ 2r and therefore ensures 
that W- dies away at a rate inversely proportional to the background shear anisotropy. This 
in turn implies that W- oc t~ 2( - 2+r \ since 9 oc i -1 . The linear decay of the W_ parameter is in 
agreement with the stability of the vacuum, plane-wave equilibrium points that consist the late 
time asymptotes of the Bianchi Vllh spacetimes (e.g. see [1, 9, 20]). 

The situation changes in the presence of matter. The latter, even when it is in the form 
of a low-density dust component, introduces new degrees of freedom into the system and the 
decrease of VV_ is not always guaranteed. Indeed, the general solution of Eq. (61) reads 

W- = dt 1 - 2 ' + C 2 t~ 2r , (63) 

where C\, C 2 are constants and C\ = Cr(l — r)(l — 2r)/(l + 2r) 4 . According to the above, 
the expansion-normalised scalar W_ will start increasing as WL oc t l ~ 2r when < r < 1/2. 
The latter corresponds to S > 1/4, which means that three quarters of the allowed Lukash 
models, those with the largest background shear anisotropy, are unstable against linear Weyl 
curvature distortions. Incidentally, the same family was also found vulnerable to linear rotational 
distortions (see section 5.2). Finally, we note that the aforementioned matter effects are sensitive 
to the precise evolution of the background model, namely to the properties of the Lukash vacuum 
solution. Given that, one should be careful before extrapolating these results to Bianchi Vllh 
models with dust. 

6 Discussion 

Bianchi models, particularly those that contain the FRW cosmologies as special subcases, are 
essential for our understanding of the large scale anisotropy of the universe. In the family of 
Bianchi spacetimes, those of type Vllh are the most general homogeneous models containing 
the spatially open FRW universe. In the absence of matter these spacetimes reduce to the plane- 
wave solution found by Doroshkevich and Lukash [2, 5]. These vacuum models also act as the 
future attractors for the non-tilted, perfect fluid Bianchi Vllh cosmologies [43]. For this reason 
the empty Vllh Lukash model has been used to study the late-time evolution of perturbed open 
FRW universes with a conventional matter content. 

Here we have engaged a mixture of covariant and orthonormal frame methods to study 
analytically the linear response of the Lukash solution to a variety of perturbations. Our results 
show that the amount of the background shear anisotropy is crucial for the stability of the 
vacuum model. More specifically, by looking into rotational or shear perturbations, we found 
that the linear instability of these distortions is more likely in models with higher background 
shear anisotropy. When dealing with vorticity perturbations, in particular, the background shear 
can force linear rotational perturbations to grow, thus mimicking the effects of a fluid with a stiff 
equation of state. Also, when the unperturbed model has the maximum allowed kinematical 
anisotropy, linear shear distortions are no longer necessarily bounded. In the latter case the 
Lukash universe perturbs away from the family of the Bianchi Vllh spacetimes. 

Our linear analysis also considered the effects of a non-zero matter component. When the 
latter was in the form of low-density dust, we found that matter had no effect on either the 
rotation or the kinematical anisotropy of the perturbed Lukash solution. However, the presence 
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of matter (even a non-relativistic pressureless fluid) plays a key role in the evolution of the Weyl 
anisotropy of the perturbed spacetime. In particular, our study showed that the introduction 
of a pressure-free component at the linear level is the catalyst that can force the vacuum Vllh 
model to diverge from its original plane-wave nature. 

In the present paper the study of the matter effects has been confined to the case of a 
pressure-free component. It is relatively straightforward to extend this formalism to include 
the effects of pressure. Generally speaking, nonzero pressure means that the matter term in 
Eq. (61) decays faster than in the case of dust. This in turn should make the matter effects 
on W_ less pronounced. However, this rather intuitive picture will be probably complicated by 
the presence of pressure gradients, which for dust are identically zero. One could also consider 
the potential implications of large-scale magnetic fields or of a non-conventional, 'dark' matter 
component. The study of magnetic fields in Bianchi VII, for example, has lead to limits on the 
strength of a possible large-scale homogeneous field more stringent than those obtained from 
standard nucleosynthesis constraints [44]. This happens because anisotropic stresses play an 
important role in the evolution of simple anisotropic universes of Bianchi type I. The anisotropic 
trace-free stress mimics the part played by the anisotropic curvature in type VII models and, in 
combination with a perfect fluid, slows the decay of the shear anisotropy in a subtle way [45]-[47]. 
The latter leads to more severe observational consequences for the CMB. A further consideration 
for future work is the role of dark energy in the universe. For a perfect fluid with n + 3p < 0, in 
violation of the strong-energy condition, the Lukash metric is unstable and approaches the flat 
Friedmann universe as t — > oo, following with the course of power-law inflation. If p = — n, then 
the dynamics approach the de Sitter universe with exponential rapidity within the event horizon 
of any geodesically moving observer [48]- [50]. This case is less interesting from a mathematical 
point of view because all distortions are rapidly inflated away. However, it is likely to be of 
considerable astronomical interest because of the growing observational evidence that inflation 
has played some role in the very early evolution of the universe and that dark energy, with 
fi + 3p < 0, is dominating the dynamics of the universe again today. Finally, we note that 
if the universe is negatively curved with compact topology there are severe constraints on the 
possibility of any homogeneous anisotropy existing in the expansion at all: all Bianchi Vllh 
universes have to be isotropic [51]- [53]. 

Acknowledgements 

The authors would like to thank Pantelis Apostolopoulos, Alan Coley, Sigbj0rn Hervik, Henk 
van Elst and John Wainwright for helpful discussions and comments. 

References 

[1] Wainwright J and Ellis G F R Eds 1997 Dynamical Systems in Cosmology (Cambridge: 
Cambridge University Press) Sections 4-9 

[2] Doroshkevich A , Lukash V N and Novikov I D 1973 Sov. Phys. JETP 37 739 

[3] Lukash VN 1974 The Evolution of Homogeneous Anisotropic Cosmological Models, Doctoral 
Dissertation (Moscow: Phys. Tech. I.) 

[4] Lukash V N 1974 Sov. Astron. 18 164 



14 



[5] Lukash V N 1976 Nuovo Cimento 35 B 268 

[6] Collins C B and Hawking S W 1973 Astrophys. J. 180 317 

[7] Barrow J D 1982 Quart. J. R. Astron. Soc. 23 344 

[8] Barrow J D and Sonoda D H 1985 Gen. Rel. Grav. 17 409 

[9] Barrow J D and Sonoda D H 1986 Phys. Rep. 139 1 

[10] Barrow J D and Tipler F J 1986 The Anthropic Cosmological Principle (Oxford: Oxford 
University Press) Section 6.11 

[11] Collins C B 1971 Commun. Math. Phys. 23 137 

[12] Bogoyavlensky O I 1985 Methods in the Qualitative Theory of Dynamical Systems in As- 
trophysics and Gas Dynamics (Springer- Verlag) 

[13] Ma P K-H and Wainwright J 1992 in Relativity Today Ed Z Perjes (Nova Science Publishers) 

[14] Coley A A 2003 Dynamical Systems and Cosmology (Dordrecht: Kluwer) 

[15] Coley A and Hervik S 2005 Class. Quantum Grav. 22 579 

[16] Wainwright J and Lim W C 2004 Preprint gr-qc/0409082 

[17] King A R and Ellis G F R 1973 Commun. Math. Phys. 31 209 

[18] Collins C B and Ellis G F R 1979 Phys. Rep. 56 65 

[19] Apostolopoulos P S 2005 Gen. Rel. Grav. to appear (Preprint gr-qc/0407040) 

[20] Hervik S, van den Hoogen R and Coley A 2005 Class. Quantum Grav. 22 607 

[21] Eardley D 1974 Commun. Math. Phys. 37 287 

[22] Ellis G F R, Uggla C and Wainwright J 1997 in Dynamical Systems in Cosmology Eds J 
Wainwright and G F R Ellis (Cambridge: Cambridge University Press) p 107 

[23] Wainwright J, Hancock M J and Uggla C 1999 Class. Quantum Grav. 16 2577 

[24] Apostolopoulos P S 2005 Class. Quantum Grav. 22 323 

[25] Hsu L and Wainwright J 1986 Class. Quantum Grav. 3 1105 

[26] Hewitt C G, Siklos S T C, Uggla C and Wainwright J 1997 in Dynamical Systems in 
Cosmology Eds J Wainwright and G F R Ellis (Cambridge: Cambridge University Press) 
p 186 

[27] Ellis G F R 1973 in Cargese Lectures in Physics Vol VI Ed E Schatzmann (New York: 
Gordon and Breach) p 1 

[28] Ellis G F R and van Elst H 1999 in Theoretical and Observational Cosmology Ed M Lachieze- 
Ray (Dordrecht: Kluwer) p 1 (Preprint gr-qc/98 12046) 



15 



[29] Penrose R 1979 in General Relativity, an Einstein Centenary Survey Eds S W Hawking 
and W Israel (Cambridge: Cambridge University Press) p 581 

Wainwright J and Anderson P J 1984 Gen. Rel. Grav. 16 609 

Gr0n and Hervik S 2001 Class. Quantum Grav. 18 601 

Barrow J D and Hervik S 2002 Class. Quantum Grav. 19 5173 

Barrow J D and Hervik S 2003 Phys. Lett. B 566 1 

Ellis G F R and MacCallum M A H 1969 Commun. Math. Phys. 12 108 

MacCallum M AH 1973 in Cargese Lectures in Physics Vol VI Ed E Schatzmann (New 
York: Gordon and Breach) p 61 

van Elst H and Uggla C 1997 Class. Quantum Grav. 14 2673 

Ellis G F R, Siklos S T C and Wainwright J 1997 in Dynamical Systems in Cosmology Eds 
J Wainwright and G F R Ellis (Cambridge: Cambridge University Press) p 11 

Barrow J D 1987 Phys. Lett. B 187 12 

Barrow J D and Cotsakis S 1988 Phys. Lett. B 214 515 

Barrow J D 1977 Mon. Not. R. Astron. Soc. 179 47P 

Barrow J D and Tsagas C G 2004 Class. Quantum Grav. 21 1773 

Hawking S W and Ellis G F R 1973 The Large Scale Structure of Space- Time (Cambridge: 
Cambridge University Press) 

Hewitt C G and Wainwright J 1993 Class. Quantum Grav. 10 99 
Barrow J D, Ferreira P G and Silk J 1997 Phys. Rev. Lett. 78 3610 
Barrow J D 1997 Phys. Rev. D 55 7451 
Barrow J D and Maartens R 1998 Phys. Rev. D 59 043502 
Barrow J D 1995 Phys. Rev. D 51 3113 
Wald R 1983 Phys. Rev. D 28 2118 

Barrow J D 1983 in The Very Early Universe Eds G Gibbons, S W Hawking and S T C 
Siklos (Cambridge: Cambridge University) p 267 

Coley A and Wainwright J 1992 Class. Quantum Grav. 9 651 

Ashtekar A and Samuel J 1991 Class. Quantum Grav. 8 2191 

Barrow J D and Kodama H 2001 Class. Quantum Grav. 18 1753 

Barrow J D and Kodama H 2001 Int. J. Mod. Phys. D 10 785 



16 



